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In this supplementary material we first provide addi- 
tional details of our integration through transients ap- 
proach to the constitutive equation and distorted mi- 
crostructure, clarifying the approximations required to 
obtain Eqs.(ll) and (12). We then present an exact result 
for the friction kernel and details of the approximations 
leading to Eqs.(13)-(15). Finally, we provide additional 
comments regarding two continuum mechanics principles 
from the standpoint of the present work. 



I. INTEGRATION THROUGH TRANSIENTS 

Central to our integration through transients approach 
is Eq.(lO) which enables averages of functions of the par- 
ticle coordinates to be calculated [7]. Applying (10) to a 
general function / of the phase space coordinates yields 
the time-dependent average 
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The time-ordered exponential function e_ generalizes 
the exponential function for (differential) operators [20| . 
Choosing / = a a p/V, where V is the system volume, 
yields directly an exact generalized Green-Kubo relation 
for the stress tensor 

<x(£) = - fdt' (TT{K(t')&}e f ]' dsnHs) &). 
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This equation is nonlinear in k(£) due to the appearance 
of f2t(£) in the exponent. The stress autocorrelation func- 
tion ( ■ ) appearing in the above expression is the overlap 
between a stress fluctuation at time £' and a second fluc- 
tuation at later time £ which has evolved under a combi- 
nation of shear-driven and internal motion. This correla- 
tion function is approximated by calculating overlaps of 
the stress fluctuations with density fluctuations describ- 
ing slow structural relaxation of the dense system. The 
lowest order non-zero projection of the dynamics is onto 
pairs of density fluctuations and is achieved using the 
pair projector 
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Insertion of this projector either side of the time-ordered 
exponential function followed by factorization of a fourth 
moment in density fluctuations into the product of second 



moments (viz. the correlators) leads to an approximate 
expression for the stress tensor. The time evolution of 
stress fluctuations is thus represented by the transient 
density correlator $k(t, £')■ Within our approach the 
static structure factor S k serves to represent the direct 
potential interactions. For glassy states ageing leads to a 
residual time dependence of Sk which cannot be resolved 
within our approach. We therefore employ an extrapola- 
tion of the equilibrium Sk when working above the glass 
transition. 

In order to calculate the distorted structure factor we 
use our general result (10) to calculate the average of 
A PkPk = PkPk - (PkPk) 
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As for the stress we can apply the projection operator 
Pi to approximate the average in the integrand, leading 
to the expression (11) for the distorted structure factor. 
We note that in calculating Sk(t; k) we require Baxter's 
result (p k p k p ) = NS (S k + ) to treat the three 
point function which arises as a result of the projection 
operator steps pll |. 



II. EQUATION OF MOTION 

In order to close our constitutive theory we require an 
equation of motion for the transient density correlator. 
The advected wavevector appearing in the definition of 
<I>k(£, £') may be reformulated as 



$k(£,£') 
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where we introduce Sil^(t) = J2i r i ' / * T (0 ' ®i- Exact 
operator identities generalizing [7] lead directly to the 
first equation of motion 
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$q(M0) + r q (i,t )*q(t,t0) 



+ / dt' M q (t, t', i )$ q (i', to) = A q (t, to). 

The initial decay rate is given explicitly by r q (i,io) = 
q 2 (t,to)/ Sg(t,t }- Exact formal expressions are obtained 
for both the generalized diffusion kernel M q (£, t',io) and 
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correction function A q (t, to). Regarding the first equa- 
tion of motion as a Volterra integral equation of the sec- 
ond kind enables the final exact form for the equation of 
motion to be formulated 
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dt' m q (t, t', t )—^(t', t ) ) = A q (i, t ). 



The generalized friction kernel and modified correction 
function are given formally by 



m q (i, t',t ) = 
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where we have introduced the operators 
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which depend upon (50 (t) = X)i r i ' ' (^» + F») 

and Q(t,to) = l-E q P q(t ,t ) ) iVS ?(t , to) <Pq(t,fa)- The irre " 
ducible part of the dynamics is contained in 



U i (t,t',t )=e f J" isnl{s ' to) . 



All results are at this stage formally exact. Approx- 
imation of the friction kernel proceeds via two steps 
and is based on the assumption that Ui(t,t' ,ta) contains 
slow dynamics only because of coupling to higher den- 
sity modes describing structural relaxation. Firstly, we 
project the average in the numerator onto density pairs 
using the time dependent projection operator 

D /. . N _ P k(t,to)^p(t,to))^£(i,t )^j3(t,to) 
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This reduces the problem to the calculation of a four 
point correlation function. The second step is to ap- 
proximate this correlator, in the spirit of quiescent mode 
coupling theory, by a product of pair correlators. The 
modified correction function makes a negligable con- 
tribution for small accumulated strains which suggests 
the approximation A(i, to) = (22J. We thus arrive at 
Eqs.(14) and (15) given in the main text. 



III. CONTINUUM MECHANICS PRINCIPLES 

The principle of material objectivity states that the re- 
lationship between the stress and strain tensors should be 
independent of the rotational state of either the sample 
or the observer [11,16,17]. This is satisfied by the Smolu- 
chowski equation which neglects inertial effects. That 
this invariance is preserved in our approximate equa- 
tions can be explicitly confirmed by considering the im- 
position of a time dependent rotation onto an arbitrary 
flow. The shear gradient, deformation gradient, left and 
right Cauchy-Green tensors in the rotating frame are thus 
given by 

k{t) = R(t)K(t)R T (t) + R(t)R T (t) 
F(t,t') = R(t)F(t,t')R T (t') 
B(t,t') = R(t)B(t,t')R T (t) 
C-\t,t') = R{t')C- x {t,t')R T {t'), 

where R(t) is a time-dependent rotation matrix. It is a 
straightforward but laborious exercise to substitute the 
transformed tensors into expressions (8), (9), (11-14) in or- 
der to obtain the required invariance result for the stress 
tensor &(t) = R(t) a(t)R T {t). 

An invariance requirement often implicitly employed in 
continuum modeling is Oldroyd's principle of local action 
[11] which states that only neighbouring particles are in- 
volved in determining the stress at any given point. The 
k integrals in (11) and (12) are insensitive to the low-/c 
behavior of the integrands. In this sense local action is 
substantiated by the present microscopic theory. 
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